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Abstract 

The concept of discrepancy plays an important role in the study of uniformity 
properties of point sets. For sets of random points, the discrepancy is a random 
variable. We apply techniques from quantum field theory to translate the problem of 
calculating the probability density of (quadratic) discrepancies into that of evaluating 
certain path integrals. Both their perturbative and non-perturbative properties are 
discussed. 
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1 Introduction 



An important actor in numerical integration is the set of points that is used. An important 
factor in the accuracy of the numerical result is the distribution of these points in the 
integration region. The relatively slow convergence with the number of points of the 
Monte Carlo method for multivariate integration has inspired a search for point sets that 
result in a faster convergence than with sets of random points. Numerical integration with 
this kind of point sets goes under the name of Quasi Monte Carlo @]. Whether a point set 
is suitable depends, of course, on the function to be integrated. Therefore, it has proven to 
be useful to assume that the integrand belongs to a certain class of functions, the problem 
class, of which it is a 'typical' member. The problem is then translated into that of an 
average- case complexity, the squared integration error made by the numerical integration 
with the given point set, averaged over the problem class 0, |^. 

If there is not enough information about the function to identify it as a member of a 
problem class, there is no choice but to look for point sets that are as uniformly distributed 
over the integration region as possible. In this analysis, measures of non-uniformity of 
point sets, called discrepancies, are needed, and have been the topic of a great number of 
publications [§,3. A class of these discrepancies, the so called quadratic discrepancies, 
appear to be identifiable as the average-case complexities mentioned before. These also 
have been the inspiration for a number of publications 0, ||, ^, |TU|, [TT]], and have been 
for this paper. 

In order to assess the quality of a certain point set and to decide whether it will do better 
than a typical 'random-point' set as in Monte Carlo, its discrepancy has to be computed 
and has to be compared with that of sets of random points. This means that one has 
to know the probability distribution of the discrepancy under sets of random points. In 

|10|, [ill the problem of calculating this distribution has been tackled for large classes of 
discrepancies. One of the techniques used in these papers was that of Feynman diagrams as 
mnemonics in the organization of the combinatorics in the calculations. This is a method 
originally introduced for calculations in quantum field theory (QFT), and besides this 
technique, other resemblances with QFT were noticed in the mentioned papers. 

In this paper we show how the problem of calculating the probability distribution of a 
quadratic discrepancy, defined as an average-case complexity, can be cast in the form of 
a problem in terms of a QFT from the start. In particular, we show how the distribution 
can be calculated as a perturbation series in on the distribution for asymptotically 
large number of points in the point set. As examples, we use the Lego problem class 
and the Wiener problem class (defined in, for example, and in this paper) to apply 
this method. 



2 General formalism 

We shall always take the integration region to be the s-dimensional unit hypercube K = 
[0, 1)^. The point set consists of N points x^, where k = 1,2, . . . ,N labels the points 
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and V = 1,2,... , s their co-ordinates. Defined as an average-case complexity on a class of 
functions : K i— > R with measure /i, the discrepancy of the point set is given by 

f 1 ^ f 

Dm = N j r/^[0]ci/i[0] , r/^vM = Jr^'Pi^k) - j <P{x)dx . (1) 

k=i -'^ 

When Xm consists of uniformly distributed random points, then the discrepancy -D at is a 
random variable with a certain probability distribution H. This probability distribution 
has been calculated for different discrepancies in various publications [|, ^, in which 
the generating function 

G{z) = E[e"^^] (2) 

has been used, where E denotes the expectation value of a random variable. This paper 
will also concentrate on the calculation of G{z). Given G, the probability density H can 
then be calculated by the Laplace transform 

+ioo 

H{DN = t) = — [ e-''G{z)dz . (3) 
2m J 

— ioo 

From now on, we assume the measure n to be Gaussian, and propose the calculation of the 
generating function from an explicit expression in terms of /i, which we will now derive. 
The integration error ?7Ar[0] can be written as a contraction r]N[(p] = J SriN{x)(l){x) dx of 
the function with a distribution given by 

1 ^ 

Svn{x) = j^J2[Six-Xk)-l] , (4) 

k=l 

where 6{x — Xk) represents the s-dimensional Dirac (5-distribution in K. In terms of the 
distribution Srjjy, the discrepancy is given by 



Dn = N 6riNix)C{x,y)5r]Niy) dxdy , (5) 
where C is the two-point Green function of the measure /i: 

C{x,y) = I ct>{x)ct>{y)d^Ji[ct>\ . (6) 



Notice that G'(O) has to be equal to one in order for the probability distribution H to 
be normalized to one. This means that / = 1 and that Eq. (|^) indeed is the proper 
definition of the two-point Green function. Because we assume the measure to be Gaussian, 
we can write for the generating function 

G{z) = I exp{zDM) dxi ■ ■ ■ dxN (7) 



^ exp ^V2zN J SriN{y)4>{y) dy^ dfj,[(j)] dxi ■ ■ ■ dx 



N 
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If now the definition of 6riN{y) is substituted, and the integrals over xi, . . . ,XAr are 
performed, we arrive at 

G{z) = j (^j e^[<^(")-^K<^fe)'^i/lrfa;^ dfi[(P] , g = ^ . (9) 

For the Lego problem class and the Wiener problem class alternative derivations are given 
in Appendix A. 

2.1 The path integral 
2.1.1 The action 

Using Eq. (P), the generating function can be written as a Euclidean path integral (cf. ||13|| ) 
with an action S given by 

S[^] = (f){x)A{x,y)(f){y) dxdy - N log (^j c^'I'^^^^-^'k^^?^)'^?^] dx^ , (10) 

where A is the symmetric linear operator with boundary conditions which is the inverse of 
the two-point Green function under the measure /x: 

A{xi,y)C{y,X2) dy = 6{xi - X2) . (11) 

K 

The two-point Green function C satisfies the boundary conditions with both of its argu- 
ments. From now on, we will assume the boundary conditions to be included in A. Formally 
this can be realized by adding linear operators with ^-distributions centered around the 
boundaries, multiplied with an arbitrary large number Q. The large numbers guarantee 
that functions which do not satisfy the boundary conditions give no contribution to the 
path integral. Notice that, because Gaussian measures are completely defined by their 
two-point Green function, A can be used as a definition of /i; functional integrals under /i 
can be written as path integrals with an action given by 



'S'oM = I I (f){x)A{x,y)(f){y)dxdy . (12) 

'K2 
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2.1.2 Gaussian measures on a countable basis 



In [jl2| it has been pointed out that a large class of quadratic discrepancies, including the 
L2-discrepancy in any dimension, can be constructed with a Gaussian measure on a class 
of functions defined by a countable set of basis functions. In this paper, we will further 



-^For example A{x,y) = Jim — ^ + MS{x — ^)^ + M6{x — ^) 5{x — y) is the linear operator — ^ 



with boundary conditions -^{k) = 4>{\) = 



M- 
d(t> I 
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only consider function classes of this kind. We assume that the members of the class can 
be written as linear combinations 

= "^(pnUnix) , 0n e R (13) 

n 

of a countable set of basis functions {«„}. Products of the basis functions are assumed to 
be integrable. In particular, we assume that the parameters 



Wn = Un{x) dx and a^.n = / Um{x)un{x) dx (14) 
exist. On such a class of functions a Gaussian measure is defined by taking 

df.[<P] = l['-^^^^^d<p^ , a„6R. (15) 

For the measure to be suitably defined, the strengths cr„ have to satisfy certain restrictions 
which can be translated into the requirement that E[D7v] exists. With this measure, the 
discrepancy becomes 

1 ^ 

Dn = —^B{Xk,Xl) , B{Xk,Xl) = '^Cr^iUniXk) - Wn){Un{xi) - Wn) ■ (16) 
k,l=l n 

A connection with the foregoing can be established with the remark that in this case, C as 
well as A can be written in terms of the basis: 

C{Xi,X2) = '^CrlUn{Xi)Un{x2) , A{xi,X2) = ^ ^ Mri(Xi)M„(x2) . (17) 
n n " 

The basis consists of the eigenfunctions of A and the strengths correspond with the eigen- 
values: 

A{x,y)un{y)dy = XnUnix) , ^ = An . (18) 

The boundary conditions are those satisfied by the basis functions. Notice that the re- 
striction to such classes of functions is equivalent with the restriction to measures defined 
with operators A that allow for a spectral decomposition in terms of their eigenfunctions 
as in Eq. ([T7|). Following the notation that is more frequently used in the path integral 
formulation of QFT, the measure can be written as 

= exp(-S'o[0])I)0 , (19) 

but whenever accurate analyses are needed, we will refer to Eq. (|T3p. 
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2.2 Perturbation theory 



The action given by Eq. ([1(]|) is highly non-local because it multiplies function values (f){x) 
and (f){y) at finite distances \x — y\. In this respect the similarity with ordinary QFT fails. 
However, because we are mainly interested in H{t) as an asymptotic expansion in 
and we assume to be very large, it can be written as a perturbation series in and 
locality is restored if the series is truncated at finite order. The zeroth-order term of the 
action will then, as in QFT, be quadratic, and the remainder, denoted by V, will have an 
expansion starting with terms of (9(0^): 

S[<P] = \ I ^{x)A,{x,y)ct){y)dxdy + V\<t)] , (20) 

with 

A,{x, y) = A{x, y) - 2z6{x -y) + 2z (21) 

and 



V[<P\ = NgJ^ + N\g^ {h - Jl) - N\og(^j^e3't> dx!^ 

= Ng''{-\j; + \j^J2-\Jl)+NO{g^) , = f <P{xY dx 

Jk 



(22) 



Perturbation theory can now be applied to calculate G{z). The total path integral (P) is 
evaluated as a perturbation series in on the norm of a Gaussian measure Hz, naively 
defined by 



dn;,[(j)\ = exp 



-i j^^ (p{x)A,ix, y)<P{y) dxdy^ Vcjy (23) 



(notice that A^ inherits the boundary conditions included in A). The series can be written 
as a diagrammatic expansion with a propagator Qz, which is the two-point function of the 
measure fiz, so 

Gz{Xl,X2) = f. (24) 

and it satisfies 

/ Az{xi,y)Qz{y,X2) dy = S{xi - X2) . (25) 
Jk 

The diagrams are a help in the organization of the terms that contribute to a given order in 
the perturbation series. These terms are proportional to moments J (j){xi) ■ ■ ■ (pixj) dfiz[4>], 
integrated over the various co-ordinates Xt, . . . ,Xj, where the order of the moment is equal 
to the order of the term in the perturbation series. 

According to the Gaussian integration rules, the zeroth order term in the series is 
proportional to (det A^)"^/^. It is, however, not clear at this point what the remaining 
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factor is and whether det Az is defined properly. To overcome this problem, we will assume 
that the operator has a spectral decomposition in terms of its eigenfunctions. Because 
Az ^ A if z ^ 0, ^Ne know that for every eigenvalue A„ of A there is an eigenvalue Xn{z) of 
Az with A„(0) = lim^^o = A„. As a result of this and the definition of /iz, the zeroth 
order term Go{z) = J(i/x^[0] is given by 

^o(-) = n J -M-lUz)<l>l)d^}j = (Uj^) ' ■ (26) 

If we denote the eigenfunction corresponding with the eigenvalue \n{z) by u^^z-i then the 
propagator is given by 

Qz{xi,X2) = ^n.^(^l)^n,z(a^2) ■ (27) 

The higher orders orders in the perturbation series consist of convolutions of the propagator, 
multiplied with the zeroth order term. Prom now on, the index z in Az and Qz will be 
omitted. 



2.3 Gauge freedom 

An interesting feature of the propagator is that it is not unique. This is a result of the fact 
that a global translation 

0c ■ (t>{x) ^ (6'c0)(x) = 4){x) + c , c e R (28) 

only results in a change of the action that is at most linear in the functions 0: 

S[0c(t>] = S[(t)]+2acx[(t>] + \a^ , (29) 

where x is such that x[6'c0] = x[0] +c and where o; > 0. In terms of the original integration 
problem this can be explained by the fact that the integration error is the same for two 
integrands that differ by a constant and the fact that // is Gaussian. As a result of this, 
and the fact that wc may assume the 'measure' Vcf) to be invariant under Oc, we can take 
any functional ^ linear in with ^[Occj)] — i[(l>] + c and write 

j exp(-5[0])P0 = j 5(^[0]-c) exp(-5[0])(icP0 (30) 

— S^m^M-SemVcP , (31) 
a J 

with 

Se[<t>] = S[(t>] - \ax^[(t>] . (32) 
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The 5-distribution in the last expression tells us that the variable ^[0] 'decouples' from its 
perpendicular directions in 0-space. This decoupling of one degree of freedom is necessary 
in order for the total integral to exist, because the new action Se is invariant under global 
translations. Now we add to this action a function F, such that 

/oo 
exp(-F(c)) dc (33) 
-oo 

exists, so that 

J exp(-5[</)])P0 = j^y^ J exp{-F{a<f)])-SemV<p . (34) 

The result is that the total integral expressed in terms of S[(j)] is, apart from some normaliza- 
tions, completely equivalent with the integral expressed in terms of F(^[0]) — ^ax^ [</>] + 5* [0]. 
Because of this freedom in the choice of the action there is a freedom in the choice of A 
and the propagator Q, which we call the gauge freedom. As usual, under different gauges, 
individual Feynman diagrams evaluate to different results, but the perturbation series as 
a whole is gauge-invariant. 

Notice that, in the orthogonal-basis picture, a change of gauge in general results in a 
change of the basis functions and the eigenvalues. However, the zeroth order term in the 
perturbation series, for example defined as in Eq. (p6[), has to be gauge invariant. 



2.4 Inst ant ons 

An expansion of the action to evaluate the generating function only makes sense when it 
is an expansion around a minimum, so that it represents a saddle point approximation of 
the path integral. Therefore, a straightforward expansion such as just proposed, which is 
in fact an expansion around the trivial solution = 0, is only correct if it is an expansion 
around the minimum of the action, that is, if the trivial solution gives the only minimum 
of the action. General extrema of the action are given by solutions of the field equation 

^A{x,y)<j){y)dy + Ng-Ng ^ e^^fa) rfy " ° ' ^^^^ 

Depending on the value of z, non-trivial solutions may also exist. At this point it can be 
said that, because A{x,y) as well as 0(x) is real, non-trivial solutions only exist if z is real 
and non-zero so that (7 G R. In the analysis of the solutions we therefore can do a scaling 
4>{x) (— > (f){x)/g so that the action for these solutions is given by 

These non-trivial solutions we call instantons (cf. [T^), although this may not be a rig- 
orously correct nomenclature, in the field theoretical sense, for all situations we will en- 
counter. Notice that instantons under different gauges only differ by a constant. This is 
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easy to see because instantons under a general gauge, characterized by F and ^, are equal 
to solutions (f)0 of the field equation obtained from the gauge invariant action Se plus a 
constant c determined by the relation F'{^[(f)0] + c) = 0. The values of z for which they 
appear and the value of the action are gauge invariant, as can be concluded from Eq. (^) 
and Eq. (||). 

If N becomes large, then the contribution of an instanton to the path integral will 
behave as e"^^''^^, where does not depend on (Notice that 0(2;) does not depend 
on N because the field equation for these rescaled functions does not depend on A^.). The 
g-7vs[(/)]_jjj^g behavior of the instanton contribution makes it invisible in the perturbative 
expansion around = 0. If is larger than zero, this will not be a problem, because 
the contribution will be very small. If, however, is equal to zero, then the contribution 
will be more substantial, and it will even explode if S[0] is negative 0. This would really 
be a major problem, if it were not for the fact that, in the cases we encounter, z has to 
be real and larger than zero for these instantons to exist, and, according to Eq. (^, we 
want to integrate G{z) along the imaginary 2;-axis. In the end, when we want to close 
the integration contour in the complex z-plane to the right, we might meet the problem 
again. However, the function we want to integrate is an expansion in of the generating 
function, which is also an expansion in z around z = that can be integrated term by 
term, and therefore we will never face the infinite instanton contributions. 

3 The Lego problem class 
3.1 Definition 

The Lego problem class is obtained by dissecting the hypercube K into M non-overlapping 
bins and by taking the characteristic functions of the bins as the basis functions. Then 
Wn is the volume of bin n and this implies that all Wn are larger than zero. The functions 
{}n and the weights w„ moreover satisfy 

M M 
^n{x){}rn{x) = Sn,m^n{x) , ^ (x) = 1 V X G K and ^ W„ = 1 . (37) 

n=l n=l 

The coefficients a„ „i are equal to WnSn,m- Notice that, for this function class, the number 
of basis functions is not only countable but even finite. 

In the following, we restrict the strengths (T„ such that all cr^Wm are equal to 1. This 
choice models functions in which the largest fiuctuations appear over the smallest intervals. 
Although not a priori attractive in many cases, this choice is quite appropriate for particle 
physics, where cross sections display precisely this kind of behavior. Moreover, the average- 
case complexity is the same as that of the x^-goodness-of-fit test. With this choice, the 

^Notice that, to be able to do make a perturbation series around = 0, the action has to be zero for 
this solution, for else the terms would all become zero or would explode for large N. 
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discrepancy becomes 



In [r^l it has been shown that, for asymptotically large A^, the probability distribution of 



Dj^ under truly random point sets approaches a Gaussian distribution whenever M ^ oo. 
3.2 The action 

Using Eq. @ and the relations (|37D , we can write down the generating function without 
facing any ambiguities, obtaining 



M 



1 



^(^) = n (^) / e^p(-^M) ^'^'^ ' (39) 

n=l 

with 

M M / M \ nr- 

m = -5^w^n0^ + iV^?5^^„0„-iVlog 5^w;„e^<^" , g = ^^ , (40) 

n=l n=l \n=l / 

where R*^ 3 = (0i, . . . , 0m) and the integration region extends over the whole of R^. 
In section |3]^ it will be shown that, if Re 2; < |wmin with Wmin = niin„w„, then the only 
extremal point of the action is a minimum at 0„ = 0, n = 1, . . . , M, so that the saddle 
point approximation boils down to a straightforward expansion in The action can be 
written in terms of a symmetric linear operator A and a potential as 

^[0] = i(0,A0) + r[0] , (41) 

with 

An,m = (1 - 2z)wJn,m + 2zWnW^ , (42) 

V[(P] = NgY,^^n<Pn + zY,^n<Pl-z(^Y.Wn^ny -N\0g(j2wne'^"^ , (43) 

where (■,■) stands for the canonical inproduct. In order for A to define the Gaussian 
measure well, its real part has to have M positive eigenvalues 0. We assumed that Re z G 
(0, Itymin] with Wmin < |, SO that Appendix B makes clear that the real part of A indeed 
has M positive eigenvalues. 

To calculate the zeroth order term in an expansion of the path integral, the deter- 
minant of A has to be calculated. The result is that 



M 



■^Summations without without exphcit hmits from 1 to M. 

'^A does not have to be a mapping R^^ i-^ R*^ ; because z may be complex it is a hnear transformation 
RA^ ^ C^. 
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which is precisely the moment generating function for the x^-distribution with M — 1 
degrees of freedom. Now perturbation theory can be apphed to calculate the rest of the 
expansion. Therefore Feynman diagrams can be used with a propagator given by the 
inverse matrix A"^ of A: 



2z 



■ Iz 



(45) 



3.3 Gauge freedom 



Because of the gauge freedom, the path integral can equally well be defined with an action 
given by — + S[(j)], where ^[0] = and where F can be any function 

with the restriction that I[F] (Eq. (P^D) exists. We can translate this into a freedom in the 
choice of the propagator A~^, for if we take F{C,) as a series in ^, then the quadratic part 
will change the propagator. If we, for example, take F{^) = |[1 + 2z{^ — 1)]^^ with e > 0, 



then 



2z 

An^rn ^ ^n,m(^) = (1 - 22)w„5„^„ + — W„Wm , (46) 



A-l ^ A-l 



1 - 2z 



S 2z 



Wn £{1 — 2z) + 2z 



(47) 



There are two interesting limits for e. The limit oi e ^ oo (the 'Feynman gauge') results 
in a diagonal propagator, which makes Feynman calculus easy. The limit of £ ^ (the 
'Landau gauge') results in a singular propagator, that is, ^4(0) is not defined. 

3.4 Inst anions 

We start this section with a repetition of the statement that non-trivial instanton solutions 
only exist if 2; G [0, 00) (section |2.4| ). In order to investigate the instantons in the Lego 
problem class, we analyze the action in terms of the variables Un = g<pn + 2z, that is, we 
consider the integral exp(— A^E[?/]) d^^y , with 

Ay] = z + ^^Y. ""-yn - log (E ^"^'") • (48) 

We are interested in the minima of S. The 'perturbative' minimum 0„ = 0, n = 1, . . . , M 
corresponds to ?/„ = 2z, n = 1, . . . , M, and general extrema of S are situated at points y 
which are solutions of the equations 

(y) = ^ _ = ^^„e^" , A; = 1,...,M. (49) 



dyk yk 2z 

If z is positive, e^'=/yk, and therefore yk, has to be positive for every k. The result is that 
the yk can take at most two values in one solution y (Fig. |l|). If they all take the same 
value, this value is 2z, and we get the perturbative solution. If they take two values, one 
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Figure 1: ?/_ and y+. 



of them, ?/+, is larger that 1 and the other, y_, is smaller than 1. With these results, and 
the fact that Eq. implies that 



WnVn = 2Z , (50) 



we see that there are no solutions but the perturbative one if 2z < Wmin, where Wmin = 
min„ Wn- 

In the next section, the other extremal points will be analyzed and it will appear that 
minima occur with T,[y] < 0. This means that, in the limit of cxd, the integral of 

exp(— A^E) is not defined; there is a 'wall' in the complex z plane along the positive real 
side of the imaginary axis, to the right of which the generating function is not defined. 
That this is not an artifact of our approach, can be seen in the expression of the generating 
function given by Eq. ( |101| ) in Appendix A. It is shown there that the generating function 



is not defined if Re 2; > J"'^^ log Wn for any one of the Wn- 

We know that, on the perturbative level, the generating function has a singularity at 
z = ^, but the instanton contributions cannot correspond with it, because they will appear 
already for Re z < |. However, in order to calculate the probability density H with the 
Laplace transform, using the perturbative expression of G{z), we can just calculate the 
contribution of the singularity at z = |, for that is the contribution to the perturbative 
expansion of H(t). 

3.4.1 The Wall 

To expose the nature of the extrema of E, we have to investigate the eigenvalues A of the 
second derivative matrix Z\ of E in the extremal points. This matrix is given by 

9^E 

^k,i{y) = (y) = ak{y)5k,i + hk{y)hi{y) , (51) 

oykoyi 
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with 



akiy) = - Vk) and bkiy) = . (52) 

To show that E becomes negative, we only use its minima, and these correspond with 
extremal points in which all eigenvalues of A are positive. According to Appendix B, we 
are therefore only interested in cases where the degeneracy of negative is one, for else 
A = Ofe would be a solution. We further are only interested in cases where there is only one 
negative a^, for if there where more, say and ak+i with < a^+i, then there would be 
a solution afc < A < a^+i < 0. So we see that the only extremal points we are interested in 
have all co-ordinates yk equal, or have one yk = y+ and the others equal to ?/_. If they are 
all equal, then they have to be equal to 2z, and for the extremal point to be a minimum 2z 
has to be smaller than 1. This is the perturbative minimum. Whether the other extremal 
points are minima depends on whether det[Z\] is positive in these points. The determinant 
can be written as 



deW)l ^ ik- (1 -l)(i^.<i^) 



(53) 



Now we notice that all extremal points can be labeled with a parameter v by defining 

pj/±(f) 

with V G (l,oo) . (54) 



We see that y± is a continuous and differentiable function of v and we have that dy±/dv = 
y±/{y± ~ !)• This parameterization induces a parameterization of 2z, and with the help 
of Eq. (|50D we see that 



d{2z) w+y+ {l-w+)y^ 
dv 2/+ - 1 - 1 

So we see that the sign of det[z^] is the same as the sign oi d{2z) / dv: if an extremal point is 
a minimum, then d{2z)/dv > 0. The minimal value that v can take to represent a solution 
is 1, which corresponds to = = 1 and 2z = 1. It is easy to see that d{2z)/dv — > — oo 
if f I 1 and w+ < |, where w+ is the value of the weight belonging to the co-ordinate with 
the value ?/+. This means that if v starts from v = 1 and increases, then it will represent 
solutions with d{2z)/dv < 0, which are local maxima. We know that, if f — >■ oo, then 
|/_ 0, cxD and 2z = w^y+ + (1 — w+)y_ — >• oo, so that d{2z)/dv has to become 

larger than at some point. The first point where 2z becomes equal to 1 again we call Vc, 
so 2z(fc) = 2z{l) = 1 (Fig. Also the function S itself can be written in terms of z{v) 
in the extremal points. Therefore we use that 

-^[w+yl + il-w+)y^__] = 4^ + 4^ (56) 
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Figure 2: S and z for instanton solutions parameterized with v. 



and that w+y^ + (1 ~ ^4-)^- = 1 if = 1, so that 



= z{v) + 



z(x) dx + 1 



Az{v) 



V — \og[2z{v)] . 



(57) 



Now the problem arises. From the previous analysis of z{v) we know that, ii 1 < v < Vc, 

(58) 



then z{v) < | so that 



l-Vc + 2 



z(x) dx < . 



Furthermore, we find that 



dv 



1 - 4^ K?/+ + (1 - ^+)y-) 



dz _ w+{l-w+){y+-y^f dz 

dv 4^2 ' y^^) 



so that also dHjdv < in v^- So there clearly is a region in [l,Wc] where dz/dv > and 
S(f) < 0. This means that in the region ^w^i^ < z < ^ there are instanton solutions 
with negative action. The situation is shown in Fig. |^ for w^i^ = 0.09. A region where 
dz/dv > and S{v) < is clearly visible in [l,fc]. 



4 The L2-discrepancy and the Wiener problem class 
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4.1 Definition 



The standard Lg-discrepancy [Q is defined as the squared integration error made by in- 
tegrating the characteristic functions ^dy of hypercubes [0, y^) x [0, y^) x ■ ■ ■ x [0, y^) C K 
with the point set Xn, averaged over y: 



D 



N 



N 



x^) . (60) 



The Wozniakowski lemma from Ref. states that it can be written as in Eq. (|I]), that 
is, as the squared integration error, averaged with respect to a variation of the Wiener 
sheet measure in which the functions are pinned down at x = (1,1,... ,1) rather than at 
a: = (0, 0. . . . , 0). The Wiener sheet measure itself is Gaussian with the two-point Green 
function given by 

s 

C{x,y) = J]min(x^2/^) . (61) 



In [|T2l it has been shown that it has a spectral representation in terms of a set of orthogonal 



functions on K. In the case s = 1, on which we shall concentrate here, these functions are 
the eigenf unctions of the linear operator ^-^ —0" in the space of functions : K i-^ R 
with boundary conditions 0(0) = 0'(1) = 0, and, in fact, the measure can be defined with 
an action given by S[(f)] = — | Jj^ 0(x)0"(x) dx, with 0(0) = 0'(1) = 0. Usually the measure 
is written in terms of an action 

Sq[(P] = I I ^'{xfdx , with 0(0) = , (62) 
2 Jk 

from which the other one can be obtained by partial integration. 
4.2 The action 

Eq. (y) states that the generating function of the Lg-discrepancy in one dimension is given 
by a path integral with an action S given by 

S[<P\ = ]^j^<P\xfdx + Ng j^ct>{x)dx-N\og(^j^e^'t'^^Ux^ ,g = ^, (63) 

where the functions satisfy the boundary condition 0(0) = 0. In section |4.6| it will be 
shown that if Re z < ^tt^, then the only extremal point of the action is a minimum at 
= 0, so that a saddle point expansion boils down to an expansion in 1/A^. 

■'Primes stand for derivatives. 
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4.3 Gauge freedom 

The boundary condition can be included into Eq. ( |63D by adding a term |M0(O)^ with 
M — >■ oo. If M is taken finite, then the problem class is not restricted anymore to functions 
with 0(0) = 0, but 0(0) gets a Gaussian distribution. For the discrepancy, however, this 
does not matter, because the extra functions that are admitted to the problem class differ 
from the original functions only by an integration constant. For finite M, the action 
transforms as S'[0] 5'[0] +Mc0(O) + |Mc^ under a global translation 0c- Therefore, the 
path integral can equally well be defined with an action given by 

F(e[0]) - iM0(O)2 + 5[0] , (64) 

where ^ can be any linear functional such that ^[6'c0] = ^[0] + c, and F can be any 
function with the restriction that I[F] (Eq. (p3|)) exists. This results in a replacement of 
the boundary condition 0(0) = by a boundary condition dictated by F(^[0]), and this is 
the gauge freedom. 

4.4 The zeroth order contribution 

To calculate the zeroth order contribution to the path integral, the eigenfunctions and 
the corresponding eigenvalues of A have to be found. We choose the gauge in which 
^[0] = Jj^ 0(x) dx and -F(^[0]) = iM^^[0], so that A with boundary conditions is given by 



(A0)(x) = -0"(x) - 2#(x) + 2^ / (t){y)dy , (65) 
6(x) dx = Q and 0'(1) = 0'(O) = . (66) 



K 

The eigenfunctions and the eigenvalues are given by 

Uk{,x) = V2 cos{kTcx) , Xk = k'^n'^ — 2z , fc = l,2, ... . (67) 
According to Eq. (^), the zeroth order contribution is given by 

and this is the well known expression for the generating function of the probability distri- 
bution of the L2-discrepancy in one dimension for asymptotically large A^. 

In a general gauge, with quadratic F, the operator A including the boundary conditions 
is given by 



(A0)(x) = - (p" (x) - 2z^{x) + 2z / (p{y)dy 

+ [6{x - 1) - 6{x)](j)'{x) + Mp{x) I (j){y)p{y) dy , (69) 



K 
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where p is a distribution with the only restriction that it integrates to one. Notice that, if 
p = 1 is taken with M ^ oo, then A in the previous gauge is obtained. By integration of 
the eigenvalue equation of A in the general gauge, it is easy to see that 

Mp(x) / (piy)p{y)dy = Ap(x) / (Piy)dy , (70) 
Jk Jk 

where A is the eigenvalue, so that the eigenvalue equation for A = becomes 

-(j)"{x) -2z(j){x) + 2z I (l){y)dy+[5{x-l)-5{x)](t)'{x) = . (71) 

The term with the 5-distributions just gives the boundary conditions 0'(O) = 0'(1) = 
and we see that solutions only exist for 2z = fc^vr^, k = 1,2, . . . . These are the values of 
2z for which X{z) = and this result is gauge invariant. According to Eq. (^6|), they are 
equal to the values of z for which l/Go{z) = 0. With the use of the factor theorem of 
Weierstrass (cf. |jl6l) we can give a first impulse to calculate Go{z) in the general gauge 
and write, with l/G'o(^)^ = JC^z), 

/(2.) = /(0)e«f - n |(l - Ji,) } = /(0)e<S-i)- . (72) 

We have checked in a number of gauges that /'(0)//(0) = — |, but we did not bother to 
prove it in the general gauge. 



4.5 The propagator 

We give the propagator Q in the first gauge of the previous section. It has to satisfy the 
equation {AQ){x,y) = 6{x — y) including the boundary conditions. If we choose the same 
gauge as in the beginning of the previous section, then the propagator has to satisfy, with 
u = V2z, 

x,y) + u^g{x,y) = -6{x-y) , (73) 



^1 

dx"^ 



^(0,Z/) = ^(l,l/) = , / g{x,y)dx = . (74) 
ax ax Jk 

The solution can be written down directly in terms of the eigenfunctions and the eigenvalues 
of A and is given by 



cos(A;7rx) cos{kny) 
k 

^ - 7r~ {cos[m(1 - |a; + yl)] + cos[m(1 - |x - y|)]} . (76) 

Zu sm u 



si-.y) = ^Y. 'TtiTr (75) 

k=l 
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4.6 Inst anions 



To start, we repeat that non-trivial instanton solutions only exist if z G [0, oo) (section 
|2.4|) . In order to investigate the instantons in the Wiener problem class, we analyze E[0] = 
S[(t)/g\/N, because this new action does not depend on A^: 

E[0] = (t)\xfdx + j (P{x)dx-\og(^j e*(^)rfx^ . (77) 

Extremal points of this action are solutions of the field equation 

that also satisfy the boundary conditions, for which we take 0(0) = 0'(1) = at this point. 
Because the action as well as the equation is invariant under global translations, solutions 
can always be chosen such that Jj^ exp(0(?/)) dy = 1, so that the equation becomes 

-— 0"(x) + l-e'^(^) = , with 0'(1) = O and [ e'^^'^Ux = 1 . (79) 

We must also have 0'(O) = 0. The problem is now reduced to that of the motion of a 
classical particle with a mass 1/\/Az in a potential 



f/(0) = (80) 

and the solution can be written implicitly as 

/ — dx 1 , , 

Vii— = , , 81 

d(l) - U{(f)) 

where the integration constant E, the energy, has to be larger than zero for solutions to 
exist. It is easy to see that the solutions are oscillatory and that, if (f){x) is a solution with 
one bending point, then also 

, , , I (j)(kx — p) r < X < p even , , , 

'^'«(^)= :L l~ ^p'l ' P = 0,l,...,fc-1 , (82) 

I 0(1 + p — kx) I < X < p odd 

is a solution for k = 2, 3, ... . These new solutions have the same energy, but a larger 
number of bending points, namely k, and the value of z increases by a factor k'^. Hence 
we can classify the solutions according to the energy and the number bending points. This 
classification in terms of the number of bending points is quite natural and this can best 
be understood by looking at the limit of — oo. Then, the equation becomes 

-0"(x) -2z0(x) + 2z / (l){y)dy = , (83) 

Jk 
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1 



Figure 3: Instanton solutions (f)k with E = 5.7 and number of bending points k = 1,2, 3. 



with 0(0) = 0'(1) = and the solutions are given by 

^i,{x) = ^Jl[l-cos{knx)] , 2z = k'^7c'^ , k = l,2,... , (84) 

so that the instantons are completely classified with the number of bending points k. If 
becomes finite, these solutions are deformed but keep the same value of k (Fig.^. For 
given k there are infinitely many solutions classified by E. 

We now concentrate on the instantons with one bending point, because the numerical 
value of the action is independent of the number of bending points. Those instantons are 
completely characterized by their energy. The values of z for which these instantons exist 
are defined as a function of E by Eq. (|8T|), which states that 



0- 

where 0_ and 0+ are the classical turning points. They are solutions of U{(j)±) = E with 
0_ < < 0+. In classical mechanics, T{E) is proportional to the period of a particle in 
the potential U (cf. [|15]). 

The function T cannot be expressed in terms of elementary functions, but a number of 
its properties can be derived, as we shall now discuss. For small E, a quadratic approxi- 
mation of the potential can be made with (j)± = ±\/2E with the result that 

limT(E) = 71V2 =^ \imz(E) = W . (86) 

The question is now whether z is increasing as a function of E. To calculate T{E) for large 
E, U{(f)) can be approximated by —1 — for < and by e'^ for > 0, so that 

T{E ^ 00) ^ 2VE + 1 + ^ log (yE + VE-1^ , (87) 



18 



so T{E) is clearly increasing for large E. To analyze T{E) for small E, we make an 
expansion in powers of E. Therefore, we write 



/2E 



v) - fi-v)] dv 



(88) 



where / is a continuous solution of the implicit equation 

g/M - f(v) - 1 



2 ' 



^9) 



with f{v) ~ f for small v. In [|T^ it is shown that it is given by the function values on the 
principal Riemann sheet of the general continuous solution and that is has an expansion 
/(^) = 12^=o^nV"' with the coefficients a„ given by 



1 and Or 



n+1 



n-l 



[n 



l)a„_i + ^ kakan+i- 



k=2 



for n > 1 



(90) 



and with the radius of convergence equal to v47r. If we substitute the power series into 
Eq. ( ^8]) and integrate term by term, we obtain the following power series for < 2tt: 



T{E) 



(91) 



n=l, n odd 

The first few terms in this expansion are 

2 



T{E) = nV2 



E 1 
1 + 12 + 4 



139 f E 



180 V 12 



571 fE 



2880 V 12 



+ 0{E^ 



(92) 



The asymptotic behavior of the coefficients a„ has also been determined in [|l^, with the 
result that, for large and integer k, 



ar 



X < 



f_2(-)^ if n = Ak 

if n = 4A; + 1 

-2(-)'= ifn = 4A; + 2 

^+2V2{-)'' ifn = 4fc + 3 



(93) 



The results are summarized in Fig.^. Depicted are the behavior for large E, the expansion 
for small E and a numerical evaluation of the integral of Eq. (|85|) . Notice the strong 
deviation of the expansion from the other curves for E > 27i, the radius of convergence. 
For this plot the first 50 terms were used. It appears that T is indeed an increasing function 
of 

We now turn to the analysis of the value of the action for an instanton. In the foregoing, 
we have shown for which positive values of z no instantons exist. Now we will show that the 
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Figure 4: T{E) computed by numerical integration, as an expansion around E = and as 
an approximation for large E. The expansion is up to and including 0{E^^). 



action indeed becomes negative for z positive and large enough. For an instanton solution 
with one bending point, the action is given by 



S{E) 



Ti{E) 



<i)"{x)dx+ [ (P{x)dx = E + 2^'^^^'^ 
Jo 



T{E) 



VE - u{4>) 



(94) 
(95) 



With the use of the same approximations for U{(j)) as in the derivation of Eq. (|87|), it is 
easy to see that, for large E, Ti{E) is bounded by 



-/T-, -,\s/2 21og-E 



4 

■3' 



log (yE+^/E^^ 



(96) 



so that S{E) clearly becomes negative for large E. 

To investigate the behavior of S{E) for small E, we use an expansion again. It can be 
obtained using Eq. (M) and the relation 



(97) 



—^(E) = — T(E)-E—(E) . 
dE^ ' 2 ^ ^ dE^ ' 
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Figure 5: S{E) computed by numerical integration and as an expansion around E = 0. 
The expansion is up to and including 0{E'^^) and its radius of convergence is 2n. The 
curve for large E is the upper bound of Eq. (^). 



A derivation of this relation is given in Appendix C. For E [ a quadratic approximation of 
the potential can be used in Eq. (^) and we find that Ti(0) = 0, so that the expansion 
of T{E) can be substituted in Eq. (pTl) and an expansion of Ti(E) can be obtained by 
integrating term by term. The expansions of T{E) and Ti{E) can then be used to find the 
expansion of S{E) using Eq. (p^). The first few terms are 



Ti{E) 
S{E) 



TC 



V2 



E 
' 2" 



E^ 



5 E^ 973 
+ 



3456 2488320 



0{E^ 



E^ E^ 89 E^ 

+ + + OiE" 

24 432 414720 ^ 



(98) 
(99) 



In Fig.|^, we plot S{E) as obtained from the series expansion, from the asymptotic behavior, 
and from numerical integration. The conclusion is that S{E) is always negative. 



5 Conclusion 

We have introduced the machinery of QFT to calculate the moment generating function G 
of the probability distribution H under sets of random points of a quadratic discrepancy 
D AT as a perturbation series on the generating function Go of the distribution Hq for 
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asymptotically large number of points A^. We used the fact that Dn can be defined as 
an average-case complexity over a function class and presented the formula for G itself 
as an average over that function class. We interpreted this formula as a Euclidean path 
integral and introduced the saddle point approximation to generate a perturbation series 
in This series can be seen as a diagrammatic expansion with a propagator Q, which 

we have shown to possess a gauge freedom. Furthermore, we have addressed the problem 
of phenomena, identified as instantons, that can spoil the saddle point approximation, and 
have indicated the situation in which they do not. However, we have also shown that the 
instantons can cause G{z) become undefined in certain regions of the complex 2;-plane. 

As examples, we have applied the introduced machinery to the Lego discrepancy and the 
L2-discrepancy in one dimension. We identified the gauge freedom, calculated the zeroth 
order term of the expansion and gave Q in certain gauges. The zeroth order terms, which 
give Gq, are in agreement with with earlier calculations. We have shown that instantons 
appear for both of the discrepancies and that they are no threat for the perturbation series, 
but cause G{z) to be undefined for asymptotically large N when the real part of z is larger 
then a certain positive value. For the L2-case this value is ^vr^, the smallest positive value 
of z at which Gq has a singularity. 

Results of perturbative calculations supported by the techniques put forward in this 
paper will be presented in 



Appendix A 



A more rigorous proof of Eq. (j^) for the Lego problem class goes as follows. According to 
Eq. (pHD, the discrepancy is given by 



where Sn = J2k=i^n{xk) counts the number of points Xk in bin n. If the points Xk are 
truly randomly distributed, the variables Sn are distributed according to a multinomial 
distribution so that the generating function can be written as 

where the summation is over all configurations {Sn} which satisfy '^n=i = N. Notice 
that E^e^^^^j > exp{zN/wn — zN) for every n, so that the generating function is not 
defined ii N ^ oo for the values of z with Re z > ^ log w^. Using Gaussian integration 
rules and the generalized binomial theorem, it is easy to see that Eq. ( p.01| ) can be written 

as 

E[e^^-] = e-^^fni|)7KM'"p(-^g^^^^^ rf^ , (102) 
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with g = ^/2zJN. By writing the A^-th power as a power of e and substituting ?/„ = 
(pn + Ng, Eq. (^) is obtained. 

For the Wiener problem class, we can show that there is a naive continuum limit which 
results in Eq. (|63|). We use the fact that the discrepancy can be defined as the naive 
continuum limit of 

M / N M \ 2 

= nY.^I[Y.Y. ¥nM , (103) 

p=l \fc=l n=l J 

U^)=9{^<X<^), Wr. = J\r.{x) dx , ^ = O^p < u) , CT^ = ^ . (104) 

D^j^^ is the discretized version of the L2-discrepancy, obtained when in Eq. (|60D the average 
over a finite number of points i/n, n = 1, M is taken, instead of the average over the 
whole of K. Notice that a whole class of 'discrete' discrepancies can be written as Eq. ( p.03| ), 
by choosing different expressions for the and the a p. Just like the Lego discrepancy, 
such a discrepancy can be written in terms of variables 5'„ that count the number of points 
Xk in bin n, and is given by 

^ M M 

~ ^ ^ RnmSnSm 2 ^ ^ TnSn ~\~ NU , 

n,m=l n=l 
M MM 

Rnm = a^K^K^ , Tn = '^RnmWm ) ^ = RnmWnWm ■ (105) 

p=l m=l n,m=l 

In the case of the L2-type discrepancy, the matrix R is given by Rnm = m.in{n,m)/M. 
The generating function is again given as the expectation value under the multinomial 
distribution. If we assume that the matrix R is invertible and positive definite, as it is for 
the Lg-type discrepancy, use the Gaussian integration rules and the generalized binomial 
theorem and do the appropriate co-ordinate transformations, we find 



M M / M . 

= - 5^ + v^5^u;„0„-iVlog Q];7i;„e^^" . (106) 



2 

n,m=l n=l \n=l 



For the L2-type discrepancy the inverse R ^ of the matrix R is easy to find and we get 

M M 
J2 KLMm = Mcj)l + MY,i<Pn-(l^n-lf , (107) 



n,m=l n=2 



SO that a naive continuum limit clearly produces Eq. (|5BD. Notice that the term M(j)1 gives 
the boundary condition 0(0) = 0; because M becomes large, functions with 0(0) ^ will 
not contribute to the path integral. 
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Appendix B: Matrices of the form Anm = cin^nm + ^^n^m 

The eigenvalues A of a real-valued matrix A are given by the zeros of the characteristic 
polynomial Pa- If A is an M x M matrix with matrix elements 

Anm = anSnm + (^bnbm , a^, &n G R , ^ = 1, . . . , M , € = ±1 , (108) 

then the characteristic polynomial Pa is given by 

M M ,2 

Pa{x) = QA{x)\[{an-x) , Qa{x) = l + • (109) 

0,^1 X 

n=l m=l 

For simplicity, we assume that the coefficients a„ are ordered such that ai < 02 < ■ ■ ■ < aM- 
If a number of dn coefficients a„ take the same value, that is, if a„ is dn-iold degenerate, 
then A = a„ is a [dn — l)-fold degenerate eigenvalue of A. The remaining eigenvalues are 
given by the zeros of the function Qa- Except of the poles at x = a„, n = 1, . . . , M, this 
function is continuous and differentiable on the whole of R. Furthermore, the sign of the 
derivative is equal to e. This means that for each zero A of Qa except one, there is an n, 
such that a„ < A < a„+m for the nearest and non-equal neighbor an+m of an- The one 
other zero is smaller than ai if e = — 1, and larger than if e = 1- This is easy to see 
because lim^^^QA{x) = lim^^_^ Qa{x) = 1. 

Appendix C: Derivation of Eq. ( [97| ) 

We use the definitions of T{E) as the r.h.s. of Eq. (^) and Ti{E) as given in Eq. (pSf): 

0+ 0+ 



Because the end points 0+ and 0_ depend on E such that E — U{(j)±) = 0, we can use 
Leibnitz's rule for differentiation under the integral sign to write 

0+ 

TiE) =2^1 - t/(0) ci0 . (Ill) 



Now we write ^/E - U{<p) = {E -e"^ + (f)+l){E -U{(f)))-^/^ and use that 1-e* = dU/d(f), 
so that 



dU , 



ET{E)+UE)- j -j=L 



- - U{<P) d<P 

\ 0- / 



(112) 



But the last integral is equal to zero, and as a result, we obtain Eq. 
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